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Particle production in Au+Au collisions has been measured in the PHOBOS experiment 
at RHIC for a range of collision energies. Three empirical observations have emerged 
from this dataset which require theoretical examination. First, there is clear evidence of 
limiting fragmentation. Namely, particle production in central Au + Au collisions, when 
expressed as dN/drj' (rj' = t] — ybeam), becomes energy independent at high energy for 
a broad region of r]' around rj' = 0. This energy- independent region grows with energy, 
allowing only a limited region (if any) of longitudinal boost-invariance. Second, there is a 
striking similarity between particle production in e"'"e~ and Au + Au collisions (scaled by 
the number of participating nucleon pairs). Both the total number of produced particles 
and the longitudinal distribution of produced particles are approximately the same in 
e"'"e~ and in scaled Au + Au. This observation was not predicted and has not been 
explained. Finally, particle production has been found to scale approximately with the 
number of participating nucleon pairs for (Npart) > 65. This scaling occurs both for the 
total multiplicity and for high particles (3 < Pj, < 4.5 GeV/c). 
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1. INTRODUCTION 
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Figure 1. Schematic diagram of the detector setup for the 2001 running period. 

The data described in this paper were taken during the Year 2000 and Year 2001 runs 
of the Relativistic Heavy Ion Colhder at Brookhaven National Laboratory. The PHOBOS 
apparatus shown schematically in Figure |I], was used to take data at three energies: 
^ys^ = 19.6, 130, and 200 GeV. A partial test detector was also used to take data at 
^ = 56 GeVi. 

The full PHOBOS apparatus comprises several subdetectors. Silicon detectors are 
employed for vertex finding, particle tracking and multiplicity measurements. This set of 
detectors has nearly full azimuthal coverage over a large pseudorapidity range \r]\ < 5.4. 
The detector setup also includes two sets of 16 scintillator counters ("paddle counters") 
located at ±3.21 m relative to the nominal interaction point along the beam (z) axis. 
These counters cover pseudorapidity in the range 3 < jr^l < 4.5 and served as the primary 
event trigger. The collision centrality is characterized by the average number of nucleon 
participants {Npart) , determined as described in Refs. . These references also tabulate 
the values and systematic errors of (Npart) for each centrality bin. 

2. SURVEY OF BASIC MEASUREMENTS 

A key goal in the design of the PHOBOS experiment was to perform a broad and 
systematic survey of hadronic particle production in heavy ion collisions. While this 
survey is still incomplete in terms of energies and species explored, significant progress has 
been made toward this goal. This section summarizes the basic measurements available 
so far: multiplicity, particle spectra and azimuthal asymmetry of particle production. 
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The pseudorapidity densities (dNch/drj) and particle yields {d?N / dydp^) given here refer 
to primary produced particles and do not include feed-down products from weak decays 
of neutral strange particles. Corrections were made for residual effects from secondary 
interactions and weak decay feed-down as well as for particles which were absorbed or 
produced in the material surrounding the collision (primarily the Beryllium beampipe 
and the magnet steel). 




Figure 2. Charged particle pseudorapidity distribution, dNch/drj, measured for Au + Au 
at = 19.6, 130, and 200 GeV for the centrality bins: 0-6%, 6-15%, 15-25%, 25- 

35%, 35-45%, 45-55%. The most peripheral bin is excluded for the lowest energy. The 
statistical error is negligible. The typical systematic error is shown as a 90% C.L. band 
for selected centrality bins. 



The unique features of the PHOBOS experiment relevant for the charged particle mul- 
tiplicity measurements are the large acceptance of the detector and the 12 m long, 1 mm 
thick. Beryllium beampipe which limits the absorption of low momentum particles and 
the production of secondaries. Figure ^ shows the charged particle pseudorapidity distri- 
butions {dNch/drj) measured at ^/s^ = 19.6, 130, and 200 GeV for a variety of centrality 
bins for pseudorapidity in the range —5.4 < rj < 5.4 Due to the large coverage in 
7], dNch/dr] is measured over almost the full range, except for a small missing fraction at 
very high \7]\. For central events, this missing fraction is estimated to be less than 2%. 
The multiplicity results are discussed further in Sections ^ and ^. 

Figure |^ shows the results of elliptic flow as a function of pseudorapidity from PHOBOS 
for 130 and 200 GeV. A unique feature of these results is the broad pseudorapidity reach 
covering almost 11 units of pseudorapidity and showing that there is no broad boost- 
invariant region near mid-rapidity. The error bars represent la statistical errors and the 
boxes represent 90% confidence level systematic errors for the 200 GeV data points. These 
results were measured using a hit-based analysis A new, track-based, analysis is also 
available both as a cross-check and to provide V2 as a function of particle |0. 

Figure § shows the unidentified charged particle yield measured by PHOBOS for a range 



4 



CM 



0.06 



0.05- 
0.04 - 
0.03- 
0.02- 
0.01 - 



I I I I I I I I I I I I I I I I I I I I 

PHOBOS Au-Au 




7 



o V 



(130) 



^(200) (Preliminary) 



2 



-4 



-2 







4 



6 



Figure 3. Elliptic flow as a function of pseudorapidity for Au-Au collisions at 
130 and 200 GeV. 



of centralities, with rapidity, ?/^, defined using the pion mass. These data cover rapidity 
in the range 0.2 < < 1.4 and transverse momenta in the range 0.2 < < 5.0 GeV/c, 
corresponding to six decades of change in the magnitude of the yield. Figure |] shows the 
invariant yields at very low p^, measured by PHOBOS 0. These particles, which make 
it through the beampipe and then range out in the active Silicon detectors in the field 
free region of the spectrometer, are mass-identified but not charge- identified. The yields 
are shown for (tt^), (A'^) and (p + p) measured at mid-rapidity (—0.1 <y < 0.4) in the 
transverse momentum ranges from 30 to 50 MeV/c for charged pions, 90 to 130 MeV/c 
for kaons and 140 to 210 MeV/c for protons and antiprotons for the 15% most central 
Au+Au collisions at y/s^ = 200 GeV. 

The unique features of PHOBOS relevant for the momentum-measured particles are the 
close proximity of the detector to the interaction region (10-80 cm), the precise vertex 
determination (0.1-0.3 mm), good segmentation (0.4-1 mm) of the silicon detectors in 
the bend direction and, again, the small amount of material between the interaction and 
the first layers of silicon. These features allow the measurement of particles with good 
momentum resolution over a broad range of transverse momenta from 0.03-5.0 GeV/c and 
beyond (when statistics allow), as well as providing the ability to reject most secondaries 
and decay products. 

The broad Pj, range, particularly the coverage to very low p^, will be important for 
constraining the role of dynamical processes such as rescattering and radial expansion in 
these collisions. The charged particle multiplicity (and elliptic flow) measurement spans 
roughly 11 units of pseudorapidity, 27r in azimuth, and a factor of 10 in beam energy, in 
a single experiment, providing strong constraints on any description of the initial state 
and subsequent dynamics. From this broad dataset, three empirical scaling rules have 
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Figure 4. Charged particle invariant 
yields in Au + Au collisions for a variety of 
centrality ranges. Non-central data have 
been scaled by arbitrary factors labeled on 
the plot for clarity. 



Figure 5. Identified particle invariant 
yields as a function of pt for central (0- 
15%) Au + Au collisions. The boxes show 
systematic uncertainties. 



emerged. Sections |^-^ describe each of them in turn. These rules, at the very least, 
provide a compact description of the data which must be respected by any model. At 
best, these empirical scaling rules may point the way to a more accurate description of 
the dynamics of these collisions. 

3. THE LIMITING CURVE FOR PARTICLE PRODUCTION 

The pseudorapidity distributions were shown in Figure H for three different collision 
energies {^/s^). In order to separate the trivial kinematic broadening of the dNch/drj 
distribution from the more interesting dynamics, the data for Au + Au collisions at dif- 
ferent energies can be viewed in the rest frame of one of the colliding nuclei. Such an 
approach led to the ansatz of "limiting fragmentation" [§], which successfully predicted 
the energy dependence of particle production away from mid-rapidity in hadron collisions, 
including pA and pp [jlO[- This ansatz states that, at high enough collision energy, both 
d'^N/dy'dp^ and the mix of particle species (and therefore also dN/dt]') reach a limiting 
value and become independent of energy in a region around y' ~ 0, where y' = y — ybeam 
and rapidity y = tanh~^ with the z axis defined as the beam (collision) axis. 
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Figure || shows the scaled, shifted pseudorapidity distributions dNch/dr]' / (Npartf^) [§]. 
The results are folded about mid-rapidity (positive and negative rj bins are averaged). 
The distributions are observed to be independent of collision energy over a substantial t]' 
range. This is consistent with and extends a similar observation made by BRAHMS p 
over a more restricted rj' range. Both the 19.6 and 130 GeV data reach 85-90% of their 
maximum value before deviating significantly (more than 5%) from the common limiting 
curve. 

The data presented here demonstrate that 
limiting fragmentation applies in the Au + 
Au system, and that the "fragmentation 
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region" in Au + Au is rather broad, cover- 
ing more than half of the available range 
of Tj' over which particles are produced. In 
particular, the fragmentation region grows 
significantly between 19.6 GeV and 130 GeV, 
extending more than two units away from 
the beam rapidity. Particle production ap- 
pears to approach a fixed limiting curve 
which extends far from the original beam 
rapidity, indicating that this limiting curve 
is an important feature of the overall inter- 
action and not simply a nuclear breakup 

effect. This result is in sharp contrast to Figure 6: Shifted pseudorapidity distribu- 
the boost-invariance scenario [|12[| which pre- tion, dNch/dt]' per participant pair, where 
diets a fixed fragmentation region and a rj' = rj — Uheami for Au + Au data at ^/s^ = 
broad central rapidity plateau that grows 19.6, 130, and 200 GeV. Systematic errors 
in extent with increasing energy. not shown. The curve is to "guide the eye" . 

4. SIMILARITY OF AA AND e+e" AT HIGH ENERGY 

In the upper part of Figure |^, total charged multiplicity data from pp, pp, e^e~, and 
AA (scaled by {Npart)/'^) are shown as a function of the appropriate ^/s for each system. 
See Ref. 




T3| for more details and 14 



for the non-PHOBOS data. The e"'"e~ data serves 
as a reference sample, describing the behavior of a simple color dipole system with a 
large a/s. The curve is a description of the e^e~ data, given by the functional form: 

Cas{s)^e^y^^^^ with the parameters A and B calculable in perturbative QCD and the 
constant parameter C determined by a fit to the e^e~ data ||15[. In order to compare 
them with e+e~, the pp and pp data were plotted at an effective energy ^Sg// = 
which approximately accounts for the "leading particle effect" |16]. Finally, central AA 
collisions, Au -|- Au from the AGS and RHIC, and Pb -|- Pb from CERN are shown. Over 
the available range of RHIC energies from 19.6 to 200 GeV, the Au -|- Au results are 
consistent with the e"'"e~ results, suggesting a universality of particle production at high 
energy. In addition, the Au + Au data approximately agrees with the scaled pp and pp 
data suggesting that the effective energy of a high energy AA collision is approximately 
just y/s^. The lower part of Figure shows these same results, but divided by the e^e~ 
fit. This figure illustrates the approach of AA from lower energies towards the high energy 
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Figure 7. Comparison of the total charged 
multiphcity versus coUision energies for 
AA, e"'"e~, pp, and pp data, as described in 
the text. In the upper panel, the curve is 
a perturbative QCD expression fit to the 
e+e" data. In the lower panel, the data 
have all been divided by the e"'"e~ fit. 



Figure 8. Pseudorapidity distribution of 
charged particles produced in central Au+ 
Au collisions at y/s^ = 200 GeV com- 
pared with the rapidity distribution along 
the thrust axis of particles produced in 
e+e" collisions at = 200 GeV. The 
Au + Au data are normalized by iVpart/2. 
Systematic errors are shown for the Au + 
Au data. 



scaling result. 

Figure H compares the pseudorapidity distribution for central Au + Au collisions to 
the closest available analog in e"'"e~ collisions: the thrust-axis-rapidity distribution, where 
the pion mass was used for all charged particles in e"'"e~ |T^. This shows that particle 



production in Au -|- Au and e"'"e collisions approximately agree in longitudinal particle 
distribution as well as in the total yield. 



5. PARTICIPANT SCALING OF PARTICLE PRODUCTION 

Participant scaling was first established in the context of the total charged particle 
multiplicity produced in pA collisions |T^. In high energy AA collisions, small deviations 
from participant scaling for mid-rapidity dN/ drj have been seen. For instance. Figure ^ 
shows the mid-rapidity dN/dr] scaled by {Npart)/'2 measured with the same apparatus for 
y^sj^ = 19.6, 130, and 200 GeV. Results of this type are sometimes interpreted in terms 

of a mixture of "hard" and "soft" components scaling as Npl'^^ and Npart respectively |]19| . 
This interpretation is somewhat ambiguous, however, since saturation models |]19|,^ de- 
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Figure 9. Scaled pseudorapidity density 
at mid-rapidity, dN/drj/ {Npart/'^') , as a 
function of centrality for 19.6, 130, and 
200 GeV. The curves correspond to pre- 



dictions from the saturation model [19,2C . 



Figure 10. Total charged multiplicity nor- 
malized by (Npart/^) as a function of Npart 
for AA. Also shown are results for e"'"e~ 
and pp data. 



scribe the data more economically, as can be seen in Figure ^, and since both dynamic 
and kinematic effects can shift particles around in pseudorapidity as centrality changes. 

Due to the large acceptance of the PHOBOS detector, we can integrate the total charged 
particle production in these colhsions. Figure |1^ shows Nch/ {Npart f^) as a function of 
(Npart) for the three different energies. Also shown are the equivalent results for pp and 
pp as well as e"'"e~ collisions at the same values of ^i. It should be noted that the pp/pp 
data are shown at their collision y/s for this plot and not y/Sgff. Three conclusions can 
be drawn from this result. First, total charged multiplicity in high energy AA collisions 
approximately scales with wounded nucleons (Npart)- Second, total particle production in 
AA collisions per participant pair is the same as the total particle production in e"'"e~ at 
the same energy for a broad range of centralities, not just for central collisions. Finally, 
particle production in pp and pp collisions is reduced compared to AA and e"'"e~. 

While the overall charged particle production scales with Npart, at high p^, particle 
production is expected to scale with Ncoii oc Npa^^. Deviations from this scaling at high 
Pj, are likely to indicate high density effects in the initial or final state. PHENIX and 
STAR have already shown that high p^, hadron production is suppressed with respect to 
Ncoii X o'pp^hx [0- Recalling that Figure indicates that pp collisions may not be the 
ideal reference sample, we investigate the scaling of high hadron production in AA 
collisions internally, using a mid-central Au + Au collision as a reference rather than pp 
data. 

Figure |ll| shows (PN/dydp^/ {Npart/2) as a function of (Npart) scaled by the value of 
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the most peripheral point for 0.5 < y < 
bins used here are specified in Figure ^. Tl: 

{Ncoul Npart) / {Ncoll/Npart)periph.- This CUrve 

Au + Au colUsions held true over the cen- 
trality range shown. This ratio varies by a 
factor of two over this centrality range, and 
yet the particle production approximately 
scales with Npart for most bins. This 
confirms the empirical observation that par- 
ticle production at high is suppressed 
with respect to the "hard scaling" expec- 
tation. However, we note the additional 
intriguing fact that, for Npart > 65, we see 
approximate Npart scaling of both low and 
high particle production. Mid-range 
particles with p^, ~ 1.5 GeV/ c show a small 
violation of Npart scaling, but are still sup- 
pressed with respect to the naive scaling 
expectations. This particular form of "high 
Pj, suppression" could be an indication of 
initial state suppression (e.g. parton satu- 
ration) or that the final state suppression 
(e.g. "jet quenching") reaches a geometric 
maximum involving one power of length 
scale Rau oc Npirt- Of course, this appar- 
ent Npart scaling could also be accidental. 
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Figure 11: Charged hadron yield in Au + Au 
per participant pair {{Npart) 1'^), normalized 
to the yield of the most peripheral bin (45- 
50%) as a function of (Npart) for different p^.. 
The grey bands represent the 90% C.L. sys- 
tematic error. The solid line in the lower 
right-hand panel shows the expectation from 
Ncoii scaling. 



6. SUMMARY 

The PHOBOS experiment has measured a systematic dataset of Au+Au collisions at 
the RHIC collider. These data span an energy range from y/s^ = 19.6 to 200 GeV, a 
pseudorapidity range from -5.4 to +5.4 and a centrality range from 65 to 340 participating 
nucleons. It also includes particle spectra covering the Pj, range from 0.03 to 5.0 GeV/c. 

Three empirical observations have emerged from this dataset. First, there is clear ev- 
idence of limiting fragmentation in Au + Au collisions: an energy-independent region 
of T]' (rj' = T] — Hbeam)- This energy- independent region grows with energy, allowing 
only a limited region (if any) of longitudinal boost-invariance. Second, there is a strik- 
ing, unexplained, similarity of particle production in e"'"e~ and particle production per 
participating nucleon pair in Au + Au collisions. Finally, particle production has been 
found to approximately scale with number of participating nucleon pairs for (Npart) > 65. 
This scaling occurs for the total multiplicity and also for high p^ particles (including 
3 < 4.5 GeV/c). 

These empirical observations serve, at least, to characterize heavy ion collisions in an 
economical way and to challenge models. At best, the observed scaling and universalities 
may point the way to a fairly simple partonic description of this high density matter and 
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thus advance our knowledge of the strong interaction. 
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